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ABSTRACT 



We revisit the question of the nature of ultraluminous X-ray sources (ULXs) through a 
detailed investigation of their spectral shape, using the highest quality X-ray data available 
in the XMM-Newton public archives (> 10, 000 counts in their EPIC spectrum). We con- 
firm that simple spectral models commonly used for the analysis and interpretation of ULXs 
(power-law continuum and multi-colour disc blackbody models) are inadequate in the face 
of such high quality data. Instead we find two near ubiquitous features in the spectrum: a 
soft excess and a roll-over in the spectrum at energies above 3 keV. We investigate a range 
of more physical models to describe these data. Slim discs which include radiation trapping 
(approximated by a p-free disc model) do not adequately fit the data, and several objects give 
unphysically high disc temperatures (fcTin > 3 keV). Instead, disc plus Comptonised corona 
models fit the data well, but the derived corona is cool, and optically thick (r ~ 5 — 30). 
This is unlike the t ^ 1 coronae seen in Galactic binaries, ruling out models where ULXs 
are powered by sub-Eddington accretion onto an intermediate mass black hole despite many 
objects having apparently cool disc temperatures. We argue that these observed disc temper- 
atures are not a good indicator of the black hole mass as the powerful, optically thick corona 
drains energy from the inner disc, and obscures it. We estimate the intrinsic (corona-less) disc 
temperature, and demonstrate that in most cases it lies in the regime of stellar mass black 
holes. These objects have spectra which range from those similar to the highest mass accre- 
tion rate states in Galactic binaries (a single peak at 2-3 keV), to those which clearly have 
two peaks, one at energies below 1 keV (from the outer, unComptonised disc) and one above 
3 keV (from the Comptonised, inner disc). However, a few ULXs have a significantly cooler 
corrected disc temperature; we suggest that these are the most extreme stellar mass black hole 
accretors, in which a massive wind completely envelopes the inner disc regions, creating a 
cool photosphere. We conclude that ULXs provide us with an observational template for the 
transition between Eddington and super-Eddington accretion flows, with the latter occupying 
a new ultraluminous accretion state. 
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1 INTRODUCTION 

Ultraluminous X-ray sources (hereafter ULXs) are bright X-ray 
sources with Lx > 10"^^ erg s~^. These objects are not asso- 
ciated with the nuclei of galaxies, so are not powered by accre- 
tion onto a central super-massive black hole, but are too bright 
for sub-Eddington accretion onto stellar-mass black holes that ra- 
diate isotropically. An obvious solution is therefore an object in- 
termediate in mass between that of stellar mass and super-massive 
black holes: intermediate mass black holes (IMBHs), of mass ~ 
lO'^ - 10'' Mq, accreting at sub-Eddington rates (e.g. Colbert & 
Mushotzky 1999; Miller & Colbert 2004). Various pathways for the 
formation of IMBHs have been suggested, both primordial (from 
the collapse of Population III stars; Madau & Rees 2001), and on- 



going (runaway mergers of massive stars, and their subsequent col- 
lapse to a massive black hole, in the densest star-formation regions 
e.g. Portegies-Zwart et al. 2004), but problems persist in producing 
sufficient numbers of IMBHs to tally with those seen in the lo- 
cal Universe (Madhusudhan et al. 2006). Thus an alternative must 
be considered: that the majority of ULXs are single stellar rem- 
nant black holes that accrete material from a close companion star, 
probably at super-Eddington rates and/or with their X-ray emis- 
sion subject to some degree of geometric beaming (e.g. King et 
al. 2001, Begelman 2002; see also Roberts 2007 and references 
therein). More extreme beaming from relativistic jets seems un- 
likely due to the number of unbeamed sources that would need to 
be present in systems such as the Cartwheel galaxy (King 2004) 

Either conclusion remains controversial, and the question of 
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what underlies ULXs can only be ultimately solved by a direct mass 
measurement based on constraints placed by the binary orbit. Un- 
til this is achieved, more indirect methods of determining the mass 
must be used in an attempt to constrain the nature of these accret- 
ing systems. A common method is to use the temperature of the 
accretion disc emission, together with its luminosity, to determine 
the mass as Shakura-Sunyaev models predict the disc temperature 
kT ~ (M/lOM0)-i/'*(L/LEdd)^^'' keV (Shakura & Sunyaev 
1973). The highest signal-to-noise spectra of ULXs can often be 
well fitted by composite models of a disc together with a hard (F < 
2) tail, with the low measured disc temperature of ~ 0.2 keV im- 
plying a high black hole mass of ~ 10^ Mq (e.g. Miller et al. 2003; 
Kaaret et al. 2003; Miller, Fabian & Miller 2004). This supports the 
IMBH interpretation, and the resulting very sub-Eddington accre- 
tion rates (< O.lLEdd) would most likely be associated with the 
low/hard state in black hole binaries (BHBs), consistent with the 
observed hard tail. 

However, it is also clear from many studies of BHBs that disc 
models only give reliable results when the X-ray spectrum is dom- 
inated by this component (the high/soft or thermal dominant state) 
(Done & Kubota 2006). The derived disc temperature is increas- 
ingly distorted as the tail increases in importance, with very low 
temperature, high luminosity discs (implying a much larger mass 
black hole than is known to be present!) seen when both tail and 
disc are strongest, as in the very high or steep power law state 
(Kubota & Done 2004, Done & Kubota 2006). The ULX spectral 
decompositions of Miller et al. (2004) all have strong tails (> 80% 
of the 0.3 - 10 keV flux; Stobbart, Roberts & Wilms 2006, hereafter 
SRW06), so a straightforward interpretation of the derived disc pa- 
rameters is unlikely to give a robust mass estimatofl 

However, the best data also shows the limitations of this sim- 
ple spectral fitting. There is clear evidence for curvature of the 
tail at the highest energies, with a deficit of photons above 5 keV 
(Roberts et al. 2005; SRW06; Miyawaki et al. 2009). Such curva- 
ture is never seen in the low/hard state at these low energies. Only 
the very high state shows curvature at such low energies, but these 
generally have spectra which are steep. The temporal variability 
also does not correspond to that expected from an IMBH in this (or 
any other) state. The simple disc and tail spectral decomposition 
shows that the tail dominates the XMM-Newton band pass, so by 
analogy with the BHBs this emission should be strongly variable. 
Yet many XMM-Newton ULX light curves show only upper limits 
on the variability power consistent with the poisson noise level for 
the data (Feng & Kaaret 2005). Thus neither the spectral nor vari- 
ability properties convincingly correspond to any of the accretion 
states known in BHBs, making it unlikely that ULXs are powered 
by sub-Eddington flows onto an IMBH (Roberts 2007). 

This conclusion is strongly supported by population studies of 
ULXs. The sheer numbers observed in starburst galaxies together 
with their short lifetimes (implied by their location in regions of 
active star formation) mean that an unrealistically large underly- 
ing population of IMBHs must be present (King 2004). Instead it 



is much more likely that the bulk of the population are the most 
extreme examples of High Mass X-ray Binaries (HMXBs) which 
somehow exceed the Eddington limit. The high mass companion 
gives a natural origin for the high mass transfer rates required 
to power the observed luminosities (Rappaport, Podsiadlowski & 
Pfahl 2005). It also explains the association of ULXs with star- 
forming regions (Fabbiano, Zezas & Murray 2001; Lira et al. 2002; 
Gao et al. 2003) and the unbroken luminosity function connecting 
ULXs to the standard X-ray binary population (Grimm, Gilfanov & 
Sunyaev 2003). Indeed, direct evidence for ULXs possessing high 
mass donor stars comes from the identification of luminous, blue 
optical stellar counterparts to several nearby ULXs (e.g. Liu et al. 
2004; Kuntz et al. 2005; Roberts, Levan & Goad 2008). 

Many ULX host galaxies have sub-solar abundances (e.g. Lee 
et al 2006). The lower opacity of subsolar material means that mas- 
sive stars lose less mass through winds during their evolution, so 
potentially can collapse to form a relatively high mass black hole 
at the end of their stellar lifetime, with Mbh ^ 80 Mq (Fryer & 
Kalogera 2001; Heger et al. 2003; Belczynski, Sadowski & Rasio 
2004; Belczynski et al. 2009). One such large stellar mass black 
hole was discovered in IC 10 X-1, a Wolf-Rayet black hole binary. 
A radial velocity curve was constructed from repeated optical ob- 
servations, which provided a mass estimate of 23-34 Mq (Prest- 
wich et al. 2007; Silverman & Filippenko 2008). Nonetheless, even 
such massive stellar remnant black holes are required to be accret- 
ing at super-Eddington rates to explain the observed luminosities 
of the brightest ULXs. Thus the accretion flows in ULXs may not 
simply be scaled up versions of those seen in BHBs, as would be 
the case for the IMBlfl model where the flows are sub-Eddington. 
Instead, observation of ULXs may allow us to probe a new regime 
of accretion physics, a new "ultraluminous state" (Roberts 2007; 
Soria 2007). 

Here we revisit the question of what the X-ray spectra of 
ULXs can tell us about the nature of their accretion flows, and how 
this constrains the nature of the accreting object, by utilising the 
best data currently available in the XMM-Newton public archives. 
We choose to use the best available data as previous studies of sam- 
ples of ULXs (e.g. Berghea et al. 2008), whilst providing interest- 
ing results, are ultimately limited in the conclusions they can draw 
by the moderate signal-to-noise of many ULX spectra. By using 
only the highest quality data from the widest band pass, highest 
sensitivity instruments available we can hope to avoid the ambigu- 
ity of previous analyses, and make definitive statements on the ac- 
cretion processes in ULXs. The paper is arranged as follows. First 
we outline the sample selection, and the data reduction processes 
(Sections|2]&[3j. Next we investigate the X-ray spectra using a va- 
riety of empirical and physical models, in order to both clarify the 
morphology of ULX spectra in the putative ultraluminous state, and 
to investigate the physical processes underlying this phenomenon 
(Section |4). Finally we discuss the implications of our results for 
the nature of ULXs (Section|5j. 



^ It should also be noted that band pass is a potential problem in comparing 
ULX results with those of BHBs. Given their high fluxes. Galactic sources 
have been most commonly studied in detail using telescopes such as RXTE, 
so are generally characterised by their behaviour in the 3-20 keV regime. 
The more distant (and hence fainter) ULXs instead require the more sen- 
sitive, high spatial resolution telescopes of Chandra and XMM-Newton for 
detailed studies, whose CCD detectors cover the 0.3-10 keV band. This dif- 
ference in typical energy range between the ULX and BHB data mean that 
any comparison between the two must be made with care. 



^ Here and throughout the paper we distinguish IMBHs as BHs with 
masses > IOOMq, as suggested by the X-ray spectral modelling of Miller 
et al. (2003) and subsequent work. Such black holes are not formed at the 
endpoint of a recent single stellar evolution, but require more exotic origins 
such as stellai' mergers in a young super stai' cluster (e.g. Portegies-Zwail; 
& McMillan 2002) or formation from primordial Population III stai's in the 
high-redshift Universe (Madau & Rees 2001). 
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Table 1. The ULX sample. 



Source 



Alternative names 



RA (J2000) Dec. (J2000) A^h " 

(10^'' cm-2) 



(Mpc) 



(1039 ergs^l) 



NGC 55 ULX^ 



M33 X-8i 

NGC 1313 X-li 
NGC 1313 X-2i 

IC342 X-li" 

NGC 2403 X-l^ 

HoIIX-ll 

M81 X-61 
HoIXX-li 

NGC 4559 X-l^ 
NGC 5204 X-l^ 
NGC 5408 X-l" 



XMMU J001528.9-3913192 
NGC 55 (fi 
Source 7* 

CXOU JO 1 335 1 .0+303937^ 
NGC 598 ULXl*^ 
Source 3^ 
1X0 7** 
Source 4^ 
IXO 8** 

NGC 1313ULX39 
Source 5^ 

CXOU J034555.7+680455" 
IXO 22^ 

PGC 13826 ULX39 
CXOU J073625.5+65354012 

Source 21^^ 
NGC 2403 X29 
IXO 31* 

PGC 23324 ULXl^ 

CXOU J081928.99+704219.412 

NGC 3031 ULX 1^5 

CXOU J095532.98+690033.412 

M81 X-9I 

NGC 3031 lO^-* 

IXO 348 

H4415 
Source 17^ 
IXO 65* 

CXOU J123551.71+275604.112 

X-7I6 

IXO 77* 

CXOU J132938.61+582505.612 
Source 23*" 

J140319.606-412259.572i^ 
Source 25 



00 15 28.9 -39 13 19.1 

01 33 50.8 +30 39 37.1 

03 18 20.0 -66 29 11.0 

03 18 22.3 -66 36 03.8 

03 45 55.5 -H68 04 54.2 

07 36 25.6 -H65 35 40.0 

08 19 29.0 -H70 42 19.3 

09 55 32.9 -1-69 00 33.3 
09 57 53.2 -1-69 03 48.3 

12 35 51.7 -1-27 56 04.1 

13 29 38.6 -1-58 25 05.7 

14 03 19.6 -41 22 59.6 



1.71 

5.58 

3.90 
3.90 

31.1 

4.17 

3.42 

4.16 
4.06 



1.78(!) 

0.70(1) 

3.70(i) 
3.70(0 

3.3(//) 

4.20(0 

4.50(0 

3.63(i70 
3.55(0 



1.49 9.70(0 
1.39 4.80(0 
5.67 4.80(rv) 



1.1 

1.0 

3.7 
4.7 

2.8 

2.4 

14.4 

2.2 
7.5 

8.0 
5.3 
3.7 



Notes: "Absoiption column values taken from Dickey & Lockman (1990) using WEBPIMMS. ''Figures shown in brackets relate to following references, from 
which the assumed distance was taken: (0 SRW06, (u) Saha et al. (2002), ((if) Liu & Di Stefano (2008), (i'v) Karachentsev et al. (2002). '^observed X-ray 
luminosity (0.3-10.0 keV) based on the DKBBFTH model (see later). Numbers shown in superscript relate to the following references for source names: 
ISRWO6, ^Stobbart et al. (2004), ^Read, Ponman & Strickland (1997), "Schlegel, BaiTett & Singh (1997), ^Grimm et al. (2005), ^Liu & Mirabel (2005), 
■'Feng & Kaaret (2005), ^Colbert & Ptak (2002), ^Liu & Bregman (2005), "'Roberts & Warwick (2000), "Roberts et al. (2004), ^^Swartz et al. (2004), 
l^Schlegel & Pannuti (2003), "Radecke (1997), l^immier & Wang (2001), ^^^Vogler, Pietsch & Bertoldi (1997), ^^Kaaret et al. (2003). 



2 SOURCE SELECTION 

Following the example of SRW06, we aim to use only the high- 
est quality data publicly available from the XMM-Newton Science 
archive in order to provide the best characterisation of the 

structure of ULX spectra. We therefore choose only the best data 
sets, ULX observations with > 10, 000 accumulated EPIC counts 
500 independent spectral bins available for fitting). This re- 
striction is imposed based on the work of SRW06, whose analy- 
sis shows that this is a reasonable threshold for statistically distin- 
guishing between physically motivated models, particularly above 
2 keV. This constraint provides a sample of 12 sources, which are 
listed in Table [T] We note that this may not be an exhaustive list, 
but that the number of sources in our sample is probably sufficient 

^ See |http : //xmm. esac . esa . int/xsa/| 



to allow global trends to become apparent. Some of this sample of 
ULXs have been observed on more than one occasion, in which 
case we select the longest individual exposure to provide the clear- 
est view of their spectrum. The selected ULXs all reside within 
nearby galaxies (^10 Mpc) due to restrictions enforced by data 
quality, and vary in foreground Galactic absorption in the range 
1.39-31.1 X lO'^'^ cm~^. Their X-ray luminosities are representa- 
tive of the full ULX range, ~ 10^'' - a few 10*" erg s'V 



3 OBSERVATIONS AND DATA REDUCTION 

Data from the longest individual observation of each source in our 
sample were downloaded from the XSA. The data sets were re- 
duced using standard tools in XMM-Newton — S AS software (ver- 
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Table 2. Observation details. 



Table 3. Size of regions used to extract the spectral data, and details of the 
resultant spectra. 



Source Obs ID Date Off axis angle" Exp'' 

(arcmins) (s) 

30410 
8650 
90200 
90200 
19750 
58470 
40800 
88300 
80400 
38300 
33700 
99300 



Notes: "Off axis angle of source in the XMM-Newton EPIC field of 
view, ''Sum of good time intervals for each observation (taken from pn 
data), calculated as per in the text. 

sion We found that background flaring was severe enough 

in three cases that periods of data were lost, resulting in multiple 
exposure data sets within the same observation. Such flaring events 
took place during the observations of Holmberg II X-1 and M81 
X-6 and caused multiple exposures in the MOS detectors only. In 
each case we find that one of the MOS exposures was heavily con- 
taminated, so for these objects we only use the other MOS data, 
together with the pn data. The third observation to be affected in 
this way is that of M33 X-8. In this case we find that the pn data is 
split into 2 exposures whilst data from the MOS detectors is split 
into 3 exposures. The first exposure from each detector contains 
no usable information. On further examination we find that the sec- 
ond MOS exposure is heavily contaminated by flaring, we therefore 
only use the second pn and third MOS exposures for the observa- 
tion of this object. To remove any remaining flaring events from 
our observations we constructed good time interval (GTI) files from 
pn data using a full-field 10-15 keV background light curve and 
count rate criteria. The exact value of the count rate criteria used 
to construct the source GTI files vary according to field (typically 
excluding count rates higher than ~1 — 1.5cts^^)to provide the 
longest exposure whilst minimising contamination. Details of these 
observations are included in Table|2l with listed exposure times in- 
corporating the GTI corrections used during the reduction of the 
data. 

The source spectra were extracted from circular apertures cen- 
tred on the individual ULXs in each detector. This was straightfor- 
ward for the majority of data, but we found that NGC 1313 X-2 
was unfortunately positioned on the chip gap of the MOSl detec- 
tor so a polygonal source region was applied to optimise our data 
extraction. Background spectra were obtained from larger circular 
regions placed near to the source. Where possible, these were po- 
sitioned on the same chip and at a similar distance from the read 
out node as the source. The only exception was NGC 4559 X-1, 
where the MOS detectors were operating in small window mode, 
therefore alternative background regions were selected on a sepa- 
rate chip in a position as close as possible to that used in the pn 
detector, whilst in the case of M8 1 X-6 no data was contained in 
the MOS 1 observation. The size of the individual source and back- 

* See |http : //xmm. esac . esa . int/sas/| 



Source 


Extraction radius 


Spectrsl bills'^ 


Rate'' 




(arcseconds) 




(ct s^^) 




Source 


Background 






NGC 55 ULX 


34 


51 


884 


2.13 


M33 X-8 


50 


75 


1252 


9.56 


N LrC lilJ A- i 


40 


60 


ioio 


1 1 n 
1.19 


NGC 1313 X-2 


40* 


60 


1600 


1.04 


1C342 X-1 


36 


54 


516 


0.68 


NGC 2403 X-1 


35 


52.5 


843 


0.48 


Ho II X-1 


52 


78 


1358 


4.93 


M81 X-6 


22 


33 


989 


0.69 


Ho IX X-1 


42 


63 


2139 


2.49 


NGC 4559 X-1 


34 


51 


593 


0.50 


NGC 5204 X-1 


40 


60 


873 


1.55 


NGC 5408 X-1 


40 


60 


990 


1.42 



Notes: "Number of spectral bins available for fitting from combined EPIC 
detectors; ''combined EPIC count rate of source during observation in the 
0.3-10 kev band. * Alternative extraction region used in MOSl because 
source positioned on edge of chip gap. 



ground extraction regions are recorded in Table [5] (in the case of 
NGC 1313 X-2 we only list the size of the circular apertures used 
in pn and MOS2). 

The best quality data (FLAG = 0) were extracted in each case 
with PATTERN 4 for pn and PATTERN ^ 12 for MOS. The 
response and ancillary response files were created automatically by 
the standard XMM-Newton tasks and spectral files were grouped to 
a minimum of 20 counts per bin, to improve statistics. The number 
of independent spectral bins after this process was completed are 
listed in Table[3] along with the combined EPIC source count rates. 



4 ULX SPECTRAL PROPERTIES 

Our aims in this work are twofold. Firstly we aim to investigate the 
basic shape of the ULX X-ray spectra, and in doing so evaluate the 
evidence for the presence of a soft excess and a power-law break 
(at energies of a few keV) in ULX spectra. As these have been sug- 
gested as the two distinguishing spectral features of a new, ultra- 
luminous accretion state (Roberts 2007) it is important to examine 
the evidence for their presence in the highest quality ULX spectra 
offered by XMM-Newton, and so determine the validity of claims 
of a new accretion state. Secondly, we will investigate the physical 
insights that a range of models can afford us on the nature of the ac- 
cretion flows in this putative state. The models we use vary from the 
simplest empirical models (power-law continua), through accretion 
disc models, to models in which we consider both an accretion disc 
and a Comptonising corona surrounding its inner regions, and the 
interplay between the two. 

All spectra are fit in XSPEC version 11.3.2 over the 0.3- 
10.0 keV energy range (unless otherwise stated). To maximize data 
quality, pn and MOS data were fit simultaneously, with the addi- 
tion of a constant multiplicative factor to compensate for calibration 
difference between the cameras. The pn constant is fixed at unity, 
whilst those for each MOS camera remain free, with the fitted val- 
ues generally agreeing to within ~10 % (larger discrepancies only 
occurred with disparate extraction regions; see Section[3]above for 
more details). In each case, the spectra are fit with two absorption 



NGC 55 ULX 


0028740201 


2001- 


■11- 


■14 


4.2 


M33 X-8 


0102640101 


2000- 


■08- 


■04 


0.4 


NGC 1313 X-1 


0405090101 


2006- 


■10- 


■15 


0.8 


NGC 1313 X-2 


0405090101 


2006- 


■10- 


■15 


6.2 


IC 342 X-1 


0206890201 


2004- 


■08- 


■17 


3.6 


NGC 2403 X-1 


0164560901 


2004- 


■09- 


■12 


5.4 


Ho II X-1 


0200470101 


2004- 


■04- 


■15 


0.3 


M81 X-6 


0111800101 


2001- 


■04- 


■22 


3.2 


Ho IX X-1 


0200980101 


2004- 


■09- 


■26 


0.3 


NGC 4559 X-1 


0152170501 


2003- 


■05- 


■27 


0.3 


NGC 5204 X-1 


0405690201 


2006- 


■11- 


■19 


0.2 


NGC 5408 X-1 


0302900101 


2006- 


■01- 


■13 


0.2 
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Table 4. Simple power-law continuum and MCD spectral fits. 



0.5 1 2 

Energy (keV) 



10 



Figure 1. XMM-Newton EPIC pn data from IC 342 X-1, fit with an absorbed 
multi-colour disc model (DISKBB). We plot only the pn data, shown in 
black, which we rebin to a minimum of lOcr statistical significance, or 10 
channels per data point for clarity. Although this is statistically one of the 
better fits to a MCD model within our sample, a visual inspection quickly 
reveals the residuals (particularly evident at high energies) that explain why 
it is still rejected at high significance. 



components; one fixed at the column observed along the line of 
sight within our own Galaxy as listed in Table [T] (from Dickey & 
Lockman 1990), and a second component that is allowed to vary to 
represent any absorption within the host galaxy and/or intrinsic to 
the ULX. The absorption columns are modeled using the TBABS 
model (Wilms, Allen & McCray 2000). All quoted errors are the 
90% confidence interval for one interesting parameter. 



4.1 Single component phenomenological models 

Simple, single component models have been used rather success- 
fully to describe the featureless spectra of ULXs in the low-to- 
moderate quality data regime for some years (e.g. Humphrey et al. 
2003; Swartz et al. 2004; Feng & Kaaret 2005; Winter, Mushotzky 
& Reynolds 2006). Here we use only the highest quality data to 
perform our analysis, yet we still find the the X-ray spectra of 
these objects to appear relatively smooth and featureless. There- 
fore, the application of these simple, single continuum models is 
a good starting point in their analysis. We start by applying an 
absorbed power-law continuum (PO in XSPEC syntax) and an ab- 
sorbed multi-coloured disc blackbody (hereafter MCD) component 
(DISKBB in XSPEC ; Mitsuda et al. 1984) separately to the data. 
The resultant fits can be seen in Table|4] 

Table |4] shows that a single power-law is not a particularly 
good fit to ULX data of this high quality. Although a rough repre- 
sentation of the observed spectra (reduced < 2) is found for 
11/12 objects, only two objects provide statistically acceptable fits 
(null hypothesis probability > 5%). Notably, these two objects - IC 
342 X-I and NGC 4559 X-I - also have the worst quality data in the 
sample. The vast majority of the other ULX data reject this model 
at very high significance. The situation is even worse when fitting 
the MCD model, where the quality of the fits are so bad they can 



Source 



Model parameters 



k-2/DoFt 



TBABS*TBABS*PO 



NGC 55 ULX 


0.455±0.008 


3.30±0.03 


1288.8/879 


M33 X-8 


0.249* 


2.218* 


2527.3/1247 


NGC 1313 X-1 


0.188±0.005 


1.85±0.02 


2138.5/1611 


NGC 1313 X-2 


0.341±0.008 


I.81±0.02 


1982.4/1595 


IC 342 X-1 


0.57±0.04 


1.83±0.05 


534.3/511 


NGC 2403 X-1 


0.51±0.02 


2.38±0.03 


1247.3/838 


Ho II X-1 


0.158±0.003 


2.63±0.01 


1602.6/1363 


M81 X-6 


0.39±0.01 


2.09±0.02 


1825.9/985 


Ho IX X-1 


0.105±0.003 


l.DUD_g QQg 


2863.4/2112 


NGC 4559 X-1 


0.120±0.009 


2.29±0.04 


586.5/588 


NGC 5204 X-1 


0.153 ±0.006 


2.52±0.03 


986.3/868 


NGC 5408 X-1 


091+0-004 


3.12±0.02 


1801.9/985 



TBABS*TBABS 








*DISKBB 








NGC 55 ULX 


0.123* 


0.573* 


2067.2/879 


M33 X-8 


0.007±0.003 


1.11 ±0.01 


1470.7/1247 


NGC 1313 X-1 


0.008* 


1.420* 


5531.5/1611 


NGC 1313 X-2 


0.115±0.004 


1.53±0.02 


2091.5/1595 


IC 342 X-1 


0.19±0.02 


, 7<:+0.07 
^■'■^-0.06 


775.1/511 


NGC 2403 X-1 




1.06±0.02 


891.0/838 


Ho II X-1 


0.0* 


0.580* 


8242.0/1363 


M81 X-6 


0.111±0.006 


1.31±0.02 


1210.1/985 


Ho IX X-1 


0.0* 


1.624* 


9164.3/2112 


NGC 4559 X-1 


0.0* 


0.708* 


1545.6/588 


NGC 5204 X-1 


0.0* 


0.618* 


2533.2/868 


NGC 5408 X-1 


0.0* 


0.314* 


8294.8/985 



Notes: Model is abbreviated to XSPEC syntax: TBABS - absorption compo- 
nents for both Galactic and external absorption; PO - power-law; DISKBB 
- MCD. Specific notes: "External absorption column (x 10^-^ cm~^) left 
free during fitting. Galactic columns listed in Table [T] ^ power-law photon 
index; inner-disc temperature (keV). *Best fitting models to this data give 
a reduced greater than 2, hence we do not place constraints due to the 
paucity of the fit. ^Here and elsewhere we use 'DoF' to abbreviate the num- 
ber of degrees of freedom available when fitting a model. 



only roughly represent five of the twelve spectra, and only one data 
set (NGC 2403 X-1) does not reject this model at high significance. 

Nonetheless, it is instructive to compare the parameters de- 
rived from these fits, since lower quality data would not be able 
to show the inadequacy of these models. For the spectra that are 
best represented by an absorbed MCD (again using the x^ < 2 
criterion) we find that fcTin ~ 1.06 - 1.7 keV. These are close 
to those seen in high mass accretion rate Galactic sources in the 
high/soft (thermal dominated) state (e.g. McClintock & Remillard 
2006), which if taken at face value would indicate that the ULXs 
contain standard stellar mass black holes accreting at fairly high 
mass accretion rates. Conversely, with a power-law model, Ta- 
ble |4] shows that the photon index measured when these objects 
are represented by an absorbed power-law range from 1.6 < F < 
3.3. While the steepest spectra would correspond to the very high 
(steep power-law) state, those with with F < 2. 1 correspond in- 
stead to the low/hard state (McClintock & Remillard 2006). The 
low/hard state is only seen at very sub-Eddington mass accretion 
rates (< O.IL/Lsdd), so a naive interpretation of this spectral de- 
composition suggests that we are observing at least some IMBHs. 

Clearly the fits from these models contradict one another as 
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Figure 2. XMM-Newton EPIC pn data from NGC 5408 X-1 (black data 
points, rebinned to a minimum of 30cr significance or 30 channels), shown 
with the best fitting multi-colour disc plus power-law model (solid line). The 
contributions of the separate model components to the overall fit are over- 
plotted as grey dashed (power-law continuum) and dashed & dotted fines 
(MCD). This clearly shows the motivation for an additional soft component 
in the spectrum. 

individual ULXs cannot be both stellar and IMBH - we note that 
the ULX data that is best represented by MCD fits generally also 
fits with power-law continua with F < 2.3. It is therefore evident 
that problems can arise when we attempt to draw physical meaning 
from simple phenomenological models applied to ULX data. In the 
case of the observations in our sample we have sufficient data qual- 
ity to demonstrate that the underlying spectrum is more complex 
than either standard disc or power-law models. As we have no rea- 
son to think our ULX sample is atypical of the class, then by exten- 
sion this should apply to all ULXs, if sufficient data were available 
for them. This should render any physical conclusions drawn from 
simple model fits to low-to-moderate quality ULX data as suspi- 
cious. 

Single component phenomenological models do not provide 
enough flexibility in fitting to accurately constrain the majority of 
the data, nor do they consistently provide a physically realistic basis 
on which to interpret the data. As our next step, we therefore turn 
to more complex phenomenological models that have been used to 
characterise the spectra of BHBs and ULXs. 



4.2 Combined phenomenological models 

Here we combine the two simplest continuum models to further 
characterise the spectral shape and features of these sources. This 
combination of a disc component plus a power-law has been used 
extremely effectively in BHBs (McClintock & Remillard 2006 and 
references therein), since although these models are not terribly 
physical, they provide a good approximation to a disc with an op- 
tically thin Comptonising corona over the 3-20 keV range. It has 
also, of course, been used as the basis for claims of IMBHs under- 
lying ULXs (see introduction). Each source spectrum is initially fit 



by an absorbed power-law component, then a multi-coloured disc 
component is added to this to look for any improvement of fit. The 
results of these fits can be seen in Table|5] 

Figure |2] illustrates the need for a second spectral compo- 
nent in ULX spectra. Here we see the spectrum of NGC 5408 X- 1 
deconvolved with the best fitting absorbed MCD plus power-law 
model. It is clear from this data that a soft component is present in 
the spectrum, as an excess above a harder continuum. In fact, we 
find that all our sources show a significant improvement in with 
the addition of a disc component (Table [5j. This implies that some 
form of soft excess is ubiquitous in ULX spectra at this level of data 
quality. 

Interestingly, in only 7/12 cases is the resulting disc 'cool' i.e. 
with fcTin < 0.4 keV as in Miller et al. (2004). Instead we find 
that the temperature from the disc component appears to vary over 
a much wider range, from 0. 17 < fcTin < 1.7 keV, covering both 
the cool and standard disc temperature range. However, where the 
disc is hot - and so forms the predominant component at higher 
energies - we find that the power-law component contributes the 
soft excess, agreeing with previous work (e.g. Stobbart, Roberts & 
Warwick 2004; Foschini et al. 2004; Roberts et al. 2005). As with 
these analyses, we note that the power-law extends to lower ener- 
gies than its putative seed photons. This is not physically realistic 
and so requires more plausible physical modelling. 

Inspection of the residuals from these fits also reveals a possi- 
ble break or turn over above ~ 2 keV. In some cases this is also 
clearly visible in their spectra. We can see in Figure |2] that the 
power-law cannot model such a feature and so is plotted in an aver- 
age position whilst the data curves around it. Such a feature is not 
present in standard BHB states and so deserves further investiga- 
tion, which is carried out in the next section. 

Another point worthy of note is that the object with the highest 
inner-disc temperature is NGC 1313 X-2, a ULX that has exhib- 
ited much cooler disc temperatures in previous observations (e.g. 
0.16 keV, Miller et al. 2003) leading to the suggestion it harbours 
an IMBH. We note that a local minimum in was observed at 
~ 0.16 keV in our data, but that the global minimum in our best 
fit occurs at much higher temperatures (fcTin = 1.7 keV). It is 
clear that these results are contradictory ~ NGC 1313 X-2 cannot 
change from an IMBH to a stellar mass BHB between observations 
- demonstrating yet again that drawing physical conclusions from 
this soft component can be hazardous and should therefore be ap- 
proached with care. 

4.3 A high energy break 

We investigate the prevalance of a high energy downturn following 
the analysis of SRW06, i.e. by comparing a power-law and a broken 
power-law description of the data above 2 keV. We do not include 
absorption in these fits as it is not easily constrained by data above 
2 keV, and besides the wider band fits infer columns that have little 
effect above 2 keV. The resulting fits are shown in Table[6l 

The broken power-law is statistically preferred (> 98 % sig- 
nificance improvement in fit according to the F-test) in eleven out 
of the twelve ULXs in our sample, with break energies in the 
~ 3.5 — 7 keV range, and a typical steepening of the power-law 
slope by AF ~ 1 — 2. This near ubiquity is made all the more re- 
markable when it is considered that the one ULX without evidence 
for a break here - NGC 5204 X-1 - has displayed evidence for such 
a feature in previous observations (Roberts et al. 2005, SRW06). 
While this break is clearly expected from the 'hot disc' fits, where 
the MCD component dominates at high energies, the break is also 
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Table 5. Combined power-law plus MCD spectral fits. 



Source 



TBABS*TBABS*(P0-I-D1SKBB) Ax^ 



NGC 55 ULX 
M33 X-8 
NGC 1313 X-1 
NGC 1313 X-2 
IC 342 X-1 
NGC 2403 X-1 
Ho II X-1 
M81 X6 
Ho IX X-1 
NGC 4559 X-1 
NGC 5204 X-1 
NGC 5408 X-1 



0.46±0.02 
0.09±0.02 



0.26 
0.29 



-1-0.02 
-0.01 
-1-0.03 
-0.02 



n+0.2 
-0.1 
f, ,cj-t-0.08 
0.3S_(, 

0.116±0.007 
0.135±0.007 

"■^°-0.02 

0.09±0.01 



3.7±0.1 
2.0±0.1 



1.70 

2.0"^ 



0.068 



-1-0.005 



+0.03 
-0.02 
1-^0.2 
-0.1 
1.7±0.1 
3+0.3 
-0.4 

2.42±0.04 

2 6+0.4 

^■°-0.3 

1.46±0.02 

^+0.07 
■0.05 
0.08 
-0.09 
2.68±0.04 



2.9"^ 



2.14: 
2.18' 



7+0.03 
-0.04 
5+0.03 
-0.04 

0.23±0.01 
1.7±0.1 

1+0.1 
-0.09 
1 19+0.04 
'■'^-0.05 

0.37±0.02 

1 42+0.04 
^■^^ -0.05 



0.7?: 

1.05^ 



0.32^ 



7+0.02 
-0.01 

0.17±0.02 
0.39±0.02 

+0.007 
-0.003 



0.27^ 



0.1 



1042.7/877 
1216.0/1245 
1796.5/1609 
1599.7/1593 
518.0/509 
853.7/836 
1503.9/1361 
1093.9/983 
2440.0/2110 
528.1/586 
925.1/866 
1320.5/983 



246.06 
1311.30 
341.88 
382.77 
16.29 
393.69 
98.71 
732.01 
423.44 
58.36 
61.16 
481.41 



Notes: models are abbreviated XSPEC syntax, as per Table|4] Specific notes: "External absorption column in units of 10^^ atoms cm ^; ''power-law photon 
index; '^inner disc temperature (keV); ''x^ improvement over the absorbed power-law fit (see Table|4j, for two extra degrees of freedom. 



Table 6. A comparison of power-law to broken power-law spectal fits in the 2-10 keV band. 



Source 



PO 



BKNPOWER 



Ax^ 



l-P(F-test)-'' 



pa 



X^/DoF 



X^/DoF 



323.1/311 


106.7 


>99 


641.6/677 


214.0 


>99 


992.9/1041 


61.3 


>99 


1011.0/1032 


184.8 


>99 


273.7/267 


7.6 


>99 


335.0/333 


147.1 


>99 


750.5/793 


22.2 


>99 


554.9/536 


349.0 


100 


1561.1/1542 


111.5 


>99 


133.8/151 


7.3 


98.1 


316.7/305 


2.8 


73.9 


432.4/415 


9.6 


98.9 



NGC 55 ULX 
M33 X-8 
NGC 1313 X-1 
NGC 1313 X-2 
IC342 X-1 
NGC 2403 X-1 
Ho II X-1 
M81 X-6 
Ho IX X-1 
NGC 4559 X-1 
NGC 5204 X-1 
NGC 5408 X-1 



3.57±0.06 
2.60±0.03 
1.06±0.02 
1.91 ±0.02 

1 5-0+0.06 

2 67+0.05 
^■"'-0.06 
9 5-c;+0.02 
^■■^°-0.03 

2.31±0.03 
1.46±0.02 

9 99+0.07 
"^■^^-0.09 
2.40±0.06 

^■°^-0.06 



429.8/314 
855.6/679 
1054.2/1043 
1195.8/1034 
281.3/269 
478.1/335 
772.7/795 
903.9/538 
1672.6/1544 
141.1/153 
319.5/307 
442.0/417 



3.1±0.1 

7+0.08 
-0.1 

r, + 0.03 



2.17^ 



1-60-0.05 

1 53+°-^ 
'■•^-'-0.07 

1.53±0.07 
2.1±0.1 
2.51 ±0.04 

, 79+0.07 
'■'^-0.1 

1.38±0.02 
2.1±0.1 

2.80±0.06 



3.9±0.3 

^■"-0.3 

6 3+°-^ 

7+0.7 

-0.2 
6 7+0.7 

4.0±0.2 
5 4+0-5 
-'■^-0.6 
4.1±0.2 

9+0.3 
-0.4 
1-1 

-0.9 
,+3 
-2 
H 

-0.9 



3.7"' 



6-2: 

4.8j 
5^ 
7.1^ 



4 0+0.5 
^■^-0.4 

3.4±0.2 

+0.3 
0.5 
+0.2 

-0.1 
H 
-0.8 

4.0±0.3 

1+0.3 

-0.2 

3.4±0.2 

2.2±0.2 
a+0.8 
-0.4 
7+2 

-0.5 
7+7 
'-2 



2.6 
2.3 



2.7 



3.1^ 



2.8 
2.7"' 



Notes: models are abbreviated to XSPEC syntax: PO - as before; BKPOWER - broken power-law model. Specific notes: "Photon index from PO model, 
''photon index before the break in BKNPOWER model, '^break energy (keV), '^photon index after the break energy, "^x^ improvement over a single power-law 
fit, f statistical probability (in per cent) that the broken power-law model provides an improvement to the fit over a single PO model, from the F-test. 



seen at high significance in the 'cool disc' objects. These are the 
sources where the IMBH model is apparently favoured, implying 
that these ULX are analogous to the low/hard state observed in 
Galactic BHBs. Yet the high energy data show a break at ~ 5 keV 
which is not seen this accretion state. The ubiquity of the high en- 
ergy break shows that the ULX spectra are not well described by a 
cool disc plus power-law as expected from the IMBH model. 

Instead, these high quality data show that the disc is either hot, 
implying a stellar mass black hole (which must be accreting at or 
above Eddington in order to produce the observed luminosity) but 
with a previously unseen soft excess (an alternative would be that 
the power-law is simple broadening the disc spectrum to explain 
a known state; see later for more detailed discussions), or the disc 
is cool but with a high energy tail which is quite unlike that seen 
in the standard spectral states of BHBs (e.g. McClintock & Remil- 
lard 2006). This new combination of observational character- 
istics - a cool disc and a broken harder component - suggests 



that these ULXs are operating in an accretion state not com- 
monly seen in the Galactic BHBs. We term this new combination 
of observational characteristics the "ultraluminous state", and it 

seems most straightforward to assume that this state accompanies 
extremely high accretion rates onto stellar remnant black holes. As 
we also have no reason to think our ULX sample atypical of the 
class as a whole, we can only presume that this spectrum may be 
endemic to the majority of ULXs. 



4.4 More physical models: slim disc 

The investigations and results described above have focused on the 
application of simple phenomenological models used to charac- 
terise the shape of the observed spectra of these sources. This work 
has revealed the presence of both a soft excess and a break above ~ 
2 keV. These studies have also shown the limitations of these mod- 
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Table 7. Spectral fits for the p-free disc model. 



Source 



TBABS*TBABS*DISKPBB 

T- ft „C 
P 



X^/DoF 



NGC 55 ULX 
M33 X-8 
NGC 1313 X-1 
NGC 1313 X-2 
IC342 X-1 
NGC 2403 X-1 
Ho II X-1 
M81 X-6 
Ho IX X-1 
NGC 4559 X-1 
NGC 5204 X-1 
NGC 5408 X-1 



0.38 ± 0.01 



0.090 



+0.009 



0.18 



-0.01 

n+0.003 

'^-O.OOS 
n 9CI+0.009 

U.ZJ1_(, g-^ 

0.57±0.04 
0.27±0.03 
0.149±0.003 
0.23'' 



3+0.01 

-0.02 

0.103±0.003 

1+0.01 

-0.009 
7+0.007 
-0.004 
0.085±0.004 



0.121] 
0.147"* 



, ,,+0.06 
'■'^-0.05 

^■^°-0.04 

7.9±0.8 

fO.09 
-0.1 



2.31" 



13 



+ 7 



2.80"* 



1.24±0.06 

n+O.l 
-0.10 

^■°"-0.06 
8 9+2 

g 2+0.6 

3.0±0.2 

^■^'-0.05 



0.408±0.005 
599+0.01 

"■-^^'-0.009 
0.521±0.002 
n s»^+0-009 

U.JOJ_Q ggg 

0.523±0.007 
0.436±0.001 

Ogj+0.02 

"■°'-o.oi 

„+0.001 
"■-'•"-0.002 
n /(fi(^+0 004 

0.466_o 
446+0.002 

u.tto_o 003 

0.392±0.002 



1019.9/878 
1242.0/1246 
2113.2/1610 
1601.2/1594 
532.1/510 
850.5/837 
1557.0/1362 
1091.7/984 
2853.2/2111 
574.7/587 
977.5/867 
1855.5/984 



Notes: Model is abbreviated to XSPEC syntax: TBABS - absorption components for both Galactic and external absorption; DISKPBB - p-free disc. Specific 
notes: "External absorption column (x 10^^ cm"^) left free during fitting, Galactic columns listed in Table[T] inner-disc temperature (keV); '^p value, where 
disc temperature scales as r^P and r is the radius. 




Figure 3. XMM-Newton EPIC pn data from NGC 1 3 1 3 X-2, displayed as per 
previous figures (rebinned to a minimum of \5cf, or 15 channels). We show 
the data fit by a power-law component based on the slope and normalisation 
of the low energy (pre-break) component of the broken power-law fit to this 
data. It is very clear from this data that some form of break or curvature is 
present above 2 keV. 



Figure 4. XMM-Newton data of NGC 1313 X-1, deconvolved with the 'p- 
free' model. The data is plotted as per the previous figure, and binned to to a 
minimum of 20(t or 20 channels for clarity. This model is used as a simpli- 
fied description of the theoretically predicted 'slim disc' model. Although 
the p value (p ~ 0.5) for this data could be considered in support of a slim 
disc, the best fit provides an unrealistic disc temperature (fcTin ~ 7.9 keV). 



els at this level of data quality, especially the MCD fits. Disc spectra 
should be more complex than this simple sum-of-blackbodies ap- 
proach. Firstly, relativistic effects broaden the spectrum, so it is less 
strongly peaked than an MCD (Cunningham 1975; Ebisawa et al. 
2001; Li et al. 2005). Secondly, the intrinsic spectrum from mate- 
rial at a given radius is not a true blackbody as there is not enough 
absorption opacity to thermalise the emission at all energies (Davis 
et al 2005). This can be approximated by a colour temperature cor- 
rection (Shimura & Takahara 1995), but this is only an approxima- 
tion and the best current models of disc spectra show weak atomic 



features in the continuum emission. Hence, even relativistic, colour 
corrected blackbody discs are not an accurate description of the 
best theoretical disc spectra over the 0.3-10 keV bandpass, though 
they are normally an excellent fit to the 3-20 keV energy range 
(Done & Davis 2008). This can be seen in real data from disc dom- 
inated BHB systems e.g. the 0.3-10 keV spectrum from LMC X-3 
is better fit by full disc models than by diskbb (Davis, Done & 
Blaes 2006). However, only a few ULX fit well even to these full 
disc models (Hui & Krolik 2008). 

These models also assume that advection of radiation is not 
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important, yet in ULXs we may be observing a super-Eddington ac- 
cretion flow. The disc becomes so opticafly thick at these extreme 
mass accretion rates that the energy released in the mid-plane of 
the disc does not have time to diffuse to the photosphere. Instead, 
the photons are carried along with the inflowing material, being ad- 
vected radially rather than being radiated vertically (Abramowicz 
et al. 1988). Such radiatively inefficient, optically thick, advection 
dominated disc solutions (termed 'slim discs') are physically very 
different to the optically thin, advection dominated flows (ADAFs) 
of Narayan & Yi (1995), as the energy is advected in photons rather 
than in protons. The additional cooling from advected photons is 
most important in the inner regions of the disc, so the expected 
luminosity from each disc radius of the slim disc is progressively 
lower than that of a standard accretion disc at the same mass accre- 
tion rate (Abramowicz et al. 1988). Thus these models predict that 
the slim disc spectra are less sharply curved than that of a standard 
disc, so these generally give a better fit to the data. Additionally, the 
slim discs may extend down to smaller radii than the classic last sta- 
ble orbit due to the non-negligible pressure support, so their spectra 
can include higher temperature components than expected from a 
standard disc (Abramowicz et al. 1988, Watarai et al. 2000, but see 
also Beloborodov 1998). Such high temperature components, if in- 
terpreted as emission from a standard accretion disc, would give an 
underestimate of the black hole mass (e.g. Makishima et al. 2000). 

To take these physical differences into account we replace the 
standard MCD model with a modified version, the 'p-free' disc 
model. This allows the disc temperature to scale as r"*", where r 
is the radius and p is a free parameter Standard (MCD) discs have 
p fixed at 0.75, but increasing amounts of advection can be mod- 
eled by decreasing p, with fully advective discs having p ~ 0.5 (e.g 
Watarai et al. 2001, Vierdayanti et al. 2006, Miyawaki et al. 2006). 
Hence a p-free disc fitting with p ~ 0.5 would be indicative of 
the slim disc model, and also gives a broader spectrum which can 
approximate the effects of the relativistic effects and colour tem- 
perature correction. 

The p-free model gives a slightly better overall fit to the data 
than a single power-law, with all ULXs being at least roughly repre- 
sented with this gentle spectral curvature. Indeed, five out of twelve 
ULXs have statistically acceptable fits to this model. If we consider 
the values derived from best fits to the data we find, at first glance, 
that our results appear to be in good agreement with the slim disc 
model, with p ~ 0.4 - 0.6 (see Table |7](, but on closer inspection 
problems begin to emerge. Table |7] shows that while some of the 
derived inner disc temperatures are in the 1-3 keV range expected 
for such super-Eddington flows onto stellar mass black holes, some 
are extremely high, with four of the twelve sources providing fits 
above 6 keV (including two of the acceptable fits). Such tempera- 
tures are high even for maximally spinning black holes (Ebisawa 
et al. 2003), making these fits physically unrealistic. These high 
quality data sets also clearly show that the spectral curvature is 
more complex in many cases, even for this modified disc model. 
One clear example of this is NGC 1313 X-I (one of the unphysi- 
cally high temperature spectra) where there is a marked inflection 
present in the spectral data at ~ 2 keV which cannot be matched 
in the p-free models (Figure |4j this is also obvious in the spectra 
of several other ULXs, see Figure[8](. So, although a p-free model 
with slim disc characteristics (or indeed a full standard disc model: 
Hui & Krolik 2008) cannot be ruled out in the minority of cases, it 
does not seem to be a good explanation for ULXs as a class. The 
two component phenomenological model (MCD plus power-law) 
is clearly a better statistical description of the majority of ULX 
spectra, though the presence of a high energy break clearly also 



shows it limitations. As a next step we replace the power-law with 
a Comptonisation model, to explore the nature of the high energy 
break. 



4.5 Comptonisation Models 

We have demonstrated that the highest quality ULX data indicates 
many of these objects are in a new, ultraluminous accretion state. 
The next, obvious question is: what are the physics of the accretion 
flow producing this shape of spectrum? In order to further explore 
the nature of these systems, we now replace the power-law contin- 
uum with more physically realistic Comptonisation models. We use 
an alternative disc model, DISKPN (Gierlinski et al. 1999), which 
incorporates an approximate stress-free inner boundary condition 
as opposed to the continuous stress assumed in DISKBB. The re- 
sulting spectra differ by less than 5 % for the same temperature, 
but we use DISKPN so as to be able to directly compare our results 
with those of SRW06. 



4. 5. 1 DISKPN + COMPTT 

We initially use COMPTT (Titai'chuk 1994) to model the coro- 
nal emission. This is an analytic approximation to non-relativistic 
thermal Comptonisation which assumes that the seed photons for 
Comptonisation have a Wien spectrum. We tie the temperature of 
these seed photons to the temperature of the accretion disc (Tmax). 
In each case the redshift is fixed at zero due to the proximity of 
these systems, whilst the optical depth and plasma temperature are 
free to vary. We show the results of our spectral fitting in Table[8] 

We find in each case that there is a local minimum in space 
which corresponds to a hot, optically thin Comptonising corona 
{kTc ~ 50 keV, r < I), similar in nature to that seen in Galactic 
BHBs in classic accretion states. However, the global minimum in 
occurs for a fit that describes a cool, optically thick corona (kTe 
~ I - 3 keV and r ~ 6 - 80). This is a very different scenario to 
those of the standard black hole accretion states. We illustrate this 
in Figure [5] which shows versus r for Ho IX X-I. Here we find 
that a local minimum is observed at low optical depths, but a clear 
global minimum occurs at r ~ 9.5, an optically thick solution. 

To quantify this improvement, we have compared the fits 
achieved in each scenario and find that in ten of the twelve spectra 
the improvement in x^ is highly significant (Ax^ > 30) (see Ta- 
ble [8}, while the other two still show Ax^ > 3.9 i.e. the break is 
detected at 90% significance. To test this further we fix the temper- 
ature of the corona at 50 keV (all hot corona local minima fits lie 
at approximately this temperature within errors). We compare the 
resultant fits to those given in Table [8] using the F-test, finding that 
in ten out of twelve cases in our sample shows a > 99 % statistical 
probability of improvement when applying a cool optically thick 
corona over a hot optically thin one. The possible presence of such 
an extraordinary corona is therefore the first major clue as to the 
origin in the physical differences between sub-Eddington accretion 
states, and the ultraluminous state. 

A second characteristic found in fitting this model to our ULX 
data is that the disc temperatures are generally cool. In fact. Ta- 
ble [8] shows that 9/12 disc temperatures reside (on face value) in 
the IMBH range (Tmax < 0.5 keV), while the remaining three ob- 
jects have hotter discs that might imply stellar mass objects. Similar 
results - a combination of a cool disc, and an optically thick cool 
Comptonising corona - have been found in previous studies of indi- 
vidual sources (e.g. Ho IX X-1 by Dewangan et al. 2006; Ho II X-1 
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Table 8. Application of Comptonisation models: DISKPN+COMPTT spectral fits 



Source 



TBABS*TBABS*(DlSKPN+COMPTT) 



l-P{F-test)-f 



Tn- 



986.3/876 


252.6 


> 99 


1204.4/1244 


42.2 


> 99 


1705.7/1608 


192.0 


> 99 


1613.2/1592 


281.6 


> 99 


515.5/508 


4.4 


87.5 


844.7/835 


131.4 


> 99 


1375.4/1360 


44.2 


> 99 


1082.7/982 


56.8 


> 99 


2284.2/2109 


126.4 


> 99 


513.1/585 


32.4 


> 99 


886.2/865 


3.9 


87.5 


1240.9/982 


41.1 


> 99 



NGC 55 ULX 
M33 X-8 
NGC 1313 X-1 
NGC 1313 X-2 
IC 342 X-1 
NGC 2403 X-1 
Ho II X-1 
M81 X-6 
Ho EX X-1 
NGC 4559 X-1 
NGC 5204 X-1 
NGC 5408 X-1 



0.235±0.004 
0.041±0.003 

0.21±0.01 

Q ,95+0.004 

0.58+0-1O 

r, ,0+0.06 

"■1^-0.01 
0.033±0.002 

0.181±0.006 

0.099±0.006 

"■'-'-0.02 

QQ35+0.007 

O.029±0.005 



0.221 



+0.010 
0.006 



0.87^ 



7+0.04 
-0.2 
0.23±0.01 

Q 7 + 0.2 

"■'-0.1 

"■^"-0.08 
27+0.1 

"■■^'-0.05 

23+0.02 
"■^^-0.01 
7+0.1 
-0.2 
0.26±0.02 



0.7^ 



0.16±0.02 
0.26±0.03 

170+0.006 
"■^ '"-0.007 



0.83 



+0.05 



0.04 
+0.08 
0.03 
2.1±0.1 

1.51±0.02 

0+8 
-0.6 

0.98±0.04 

7+0.1 
-0.10 
, ,^+0.02 

'■^-'-o.oi 

7+0.1 

-0.09 
,+0.6 
'-0.3 
7+2 
-0.4 

1 5+f-3 



1.39"* 



2.78"* 



2.12"* 



2.27"* 
l.i 

2.2"^ 



9.9±0.4 



+100 

30 
+0.6 
0.5 
+0.4 
1 

+3 
7 



80 

8.5: 
15.3: 

T 

12.5±0.7 
5.5±0.2 

31^3 
9.5±0.5 

7.4i;-^ 

6"^^ 

" -3 

6 5 ^"°-** 



Notes: models are abbreviated to XSPEC syntax: DISKPN - accretion disc model, COMPTT - Comptonisation model. Specific notes: "External absolution 
column in units of 10^^ atoms cm~^, ''maximum temperature in the accretion disc (keV), '^plasma temperature in the Comptonising corona, '^optical depth 
of corona, '^■^ improvement , and statistical probability (in percent) of the fit improvement over a hot optically thin corona with fcTe = 50 keV fixed. 



Table 9. Application of Comptonisation models: DISKPN+EQPAIR spectral fits 



Source 



TBABS*TBABS*(D1SKPN+EQPAIR) 



X^/DoF 



NGC 55 ULX 
M33 X-8 
NGC 1313 X-1 
NGC 1313 X-2 
IC342 X-1 
NGC 2403 X-1 
Ho II X-1 
M81 X-6 
Ho IX X-1 
NGC 4559 X-1 
NGC 5204 X-1 
NGC 5408 X-1 



250^" "^ 

U.Z1U_Q QJ^ 

0.204±0.01 

+0.08 
-0.07 
+0.04 
-0.02 
+0.002 
-0.001 
1-0.02 
-0.03 
f0.005 
-0.004 
hO.Ol 
-0.02 

0.047±0.003 

5+0.002 
0.001 



0.64 

0.24 
0.050 

0.20 
0.109 

0.14 



0.033" 



0.253"* 



:i + 0.01 

-0.005 

0.63±0.02 



0.271 
0.37 
0.32 
0.32 

0.267 
0.30 



+0.004 
-0.010 
-0.02 
-0.04 
+0.02 
0.05 
+0.04 
0.05 
+0.007 
0.003 
0.1 
-0.02 



0.312"* 



7+0.01 
-0.006 
0.177±0.009 

0.314±0.01 



1.24±0.03 

70"'"0 03 
"■'"-0.02 

4.02±0.07 

7 To+0.05 

^■^8_Q Qg 

+0.5 
0.1 
+0.05 
0.07 
+0.009 

0.01 
+0.1 
0.6 



0.184 



+0.003 



3.2 
1.87 
1.011 

2.6 
4.99±0.08 

2.2±0.1 

, Oi:+0.07 
^•-'-'-0.05 

1 140"^°°^ 



+0.7 
2 

+0.5 
0.8 
17 q+0.5 
^ '-^-0.6 
11.1+1'' 



26.9 
13.2 



0.8 

+3 
0.3 

27.6±2 

7.00±0.02 
,+3 

-0.3 

14±2 
-1 

-2 
□+0.4 
-0.6 



11^ 



25: 

18.5"* 



10! 
12.9"* 



990.2/876 
1204.0/1244 
1711.6/1608 
1601.8/1592 
516.1/508 
849.3/835 
1391.2/1360 
1085.2/982 
2293.1/2109 
513.8/585 
894.3/865 
1244.2/982 



Notes: models are abbreviated to XSPEC syntax: DISKPN - accretion disc model, EQPAIR - Comptonisation model. Specific notes: "External absolution 
column in units of lO'^^ atoms cm'^, ''maximum temperature in the accretion disc (keV), '^ratio between the compactness of electron and the compactness of 
seed photon distribution, "^optical depth. 



in Goad et al. 2006; also SRW06). However, in interpreting the cool 
disc component we must consider the implications of the presence 
of an optically thick corona. It is likely that such a medium would 
mask the inner-most radii of the disc. It is also likely that material 
and power may be drawn from the disc to feed the corona. A com- 
bination of these factors could greatly impact the observed temper- 
ature of the accretion disc, invalidating any conclusions drawn on 
this basis. We return to this point later. 



4.5.2 DISKPN+EQPAIR 

We now apply an alternative, more physically self-consistent 
Comptonisation model to test our provisional result from the 
DISKPN-I-COMPTT model further. Here we use EQPAIR (Coppi 
1999), a model which allows both thermal and non-thermal elec- 



trons, and calculates the resulting spectrum without assuming that 
the electrons are non relativistic, and where the seed photons can 
have a disc or blackbody spectrum. We choose to use thermal elec- 
trons plus disc emission to allow for comparison to COMPTT, and 
again tie the temperature of the seed photons to that of the inner 
accretion disc. 

Comparing the fits to the data for DISKPN+COMPTT (Table 
[U with DISKPN+EQPAIR (Table[9l(, we find that similar values 
are achieved in all cases. Results once again indicate that the global 
fitting minima are achieved by a relatively cool accretion disc ( 0. 1 8 
< Tmax < 0.63 keV) with a cool optically thick corona (r ~ 7.0- 
27.6), thick enough to potentially hide the inner disc. 

Figure |6] illustrates the similarity in results between the two 
Comptonisation models. Here we plot typical values from our fits 
for both Comptonisation components in grey (COMPTT - dashed 
line; EQPAIR - dashed & dotted line) along with an average disc 



© 2009 RAS, MNRAS 000,[T]tT7] 



The Ultraluminous State 



11 




T 

Figure 5. Variation in o^^r a range of optical depths derived from tlie 
fitting of an absorbed DISKPN plus COMPTT to tlie data from Ho IX X-1. 
A local minimum is observed at lower optical depths but this clearly shows 
that the global minimum occurs at r ~ 9.5, an optically thick solution. 

shown in black. We can clearly see that each of these models ex- 
hibits similar curvature at higher energies. The difference in the 
spectrum of these coronal components lies at lower energies where 
subtle variations in the disc and absorption components can com- 
pensate for this. Each of these models can therefore explain the 
curvature observed at both higher and lower energies within our 
band-pass. The soft excess at low energies is modelled well by 
a cool accretion disc, whilst the optically thick corona causes the 
downturn at high energies (irrespective of model choice). 

When these results are compared to the findings of SRW06, 
we find that we achieve similar fits to the dat^E Their results also 
indicate the presence of a cool accretion disc (0.08 < fcTmax < 
0.29 keV) in all cases, whilst an optically thick corona is observed 
in almost all sources (r ranges from 0.2-33). It is therefore clear 
that these appear to be the physical characteristics of ULXs, and 
hence the ultraluminous state. This clearly demonstrates that we 
appear to be observing a radically different accretion flow to the 
classic accretion states viewed in Galactic BHBs. 

We should note, however, that there are at least two Galac- 
tic black hole candidates that exhibit similar (albeit less extreme) 
spectral traits when modelled similarly - an unusually cool accre- 
tion disc, and a cool optically thick corona - to the ULXs observed 
within our sample (when band-pass is considered). These objects 
are GRS 1915+105 and XTE J1550-564, and in each case these 
objects are thought to be accreting at, or above, the Eddington limit 
when displaying such traits (Kubota & Done 2004; Middleton et 
al. 2006; Ueda et al. 2009). This is further evidence to support the 
assertion that the ultraluminous state represents a super-Eddington 
accretion flow. 



4.6 Energetic disc-corona coupling 

We now move on to the final stage in our analysis to further ex- 
plore the results provided by Comptonisation models. The models 
above implicitly assume that we can still observe the disc down to 

^ We have five data sets in common with SRW06, and a further four sources 
with improved data. 




Energy (keV) 



Figure 6. Comparison of disc plus Comptonisation models used within our 
analysis. An average cool disc is plotted in black with coronal compo- 
nents representative of those found in our fits shown in grey (dashed line 
- COMPTT; dashed & dotted line - EQPAIR). Although the Comptonisa- 
tion components have different spectral sfiucture at lower energies, it is clear 
that the same cool optically thick corona is achieved at higher energies of the 
XMM-Newton band pass. The differences that are evident in the lower part 
of our energy range are accounted for by the other components in our model. 




0.5 1 2 5 10 

Energy (keV) 



Figure 7. XMM-Newton EPIC pn data from Ho IX X-1, fit with an absorbed 
ultraluminous model (DKBBFTH, in blue). We plot only the pn data, shown 
in black, which we rebin to a minimum of 15(t significance or 15 channels 
for clarity. We also plot the various components of the accretion system 
described by the of the model, once again disc components are plotted in 
red whilst the corona is in green. The visible regions of the outer disc and 
the optically thick corona are plotted with a soUd line, whilst the masked 
emission from the cooled, energetically-coupled inner disc is represented 
by a dotted line. This shows that the curvature at lower energies is due to 
emission from the outer accretion disc (we ai'e not able to directly observe 
emission from the innermost regions), and the break or turnover at higher 
energies is caused by the emission from the optically thick corona. 
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Table 10. The ultraluminous model: DKBBFTH 



Source 



TBABS*TBABS*(DKBBFTH) 



X^/DoF 







disc 


Rc/Rin" 


Rin'' 




kTj 




NGC 55 ULX 


239+? l!L. 

0.006 


38+"-;;? 

^ — 0.01 


1 5+°A 

•^—0.2 


700+?00 
" — 100 


2 0+°-? 
-0.5 


79+°'L 

—0.08 


10 6+°-3 

" -0.5 


iVUJ jtV O 


036+° °°* 




^■'-0.7 


OJ_6 


3 6+0-1 


6+30 


1.5+0.2 


NGC 1313 X-1 


0.211+°°°4 

^ — 0.005 


0.60+°°i 

" — 0.04 

i.i8+°:f 


9 O2+0.10 
— 0.07 


460+« 

— 30 


1.66+°°2 

^"—0.02 


2.19+°-^o 

^—0.09 


7.79+°'°2 

— 0.02 


NGC 1313 X-2 


21+0.01 

"■^^-0.02 


4 3+1 

^■•'-0.5 


J 30+80 


2 5+0.1 


1+6 

-'-I 


3 5+°-l 
-■-'-0.2 


IC 342 X-1 


S5+0 0S 
"■-'-'-0.03 


1.0±0.2 




,60+400 

10U_gQ 


1 6+°-2 
^■"-0.5 


2 4+7 


7 Q+0.2 
'■^-0.5 


NGC 2403 X-1 


(, 234+0.009 


0.85+°-f 


2 4+0.4 


170+40 
^ '"-30 


1 Q+0.2 


1 003+10 
i.ouj_Q 003 


□ 7+0.2 
^■'-0.4 


Ho 11 X-1 


079+0.008 
"•"'^-0.007 


30+° '°2 


5 6+1-° 


3900±100 


+0.02 
^■-"J-O.Ol 


13 2+1 
'■'■^-0.6 


1 43+0.02 
^■^-'-0.01 


M81 X-6 


0.19±0.02 


98+°°** 


3.3±0.2 


140±40 


2 15+°-2 
^■^-'-0.07 


■ ,0+0.07 
'--'^-0.2 


6 87+0.2 
"■'^'-0.07 


Ho IX X-1 


r, n, +0.008 

"•'^'-0.005 


1 01+°°3 
^■"'-0.04 


4 0+0.2 
^■"-0.1 


220+50 


1.58±0.03 


2.5±0.2 


7 02+0.03 
'■^^^-0.03 


NGC 4559 X-1 


0.138lg:°f 


3J+0.06 


2 0+0.2 

^•"-0.1 


3000+3000 

-'"""-2000 


, Q<: + 0.09 
i-'-'-O.lO 


1 Q+0.6 


6 73+0.09 
"■'-'-0.10 


NGC 5204 X-1 


0.036±0.007 


0.54±0.02 




610+100 


2.0±0.2 


I.9±0.4 


6 5+°-2 


NGC 5408 X-1 


029+0 °°'' 


n occ +0.006 
u.zjj _0.005 


1 47+0.01 
^■^'-0.05 


3000^4°° 


9 T,+0.04 
^■-' -0.03 


1 6+°-2 


S Rn+0-04 
■'■^"-0.03 



989.1/876 
1201.9/1244 
1708.5/1608 
1594.8/1592 
515.7/508 
846.9/835 
1399.4/1360 
1080.7/982 
2286.9/2109 
513.8/585 
889.1/865 
1246.5/982 



Notes: models are abbreviated to XSPEC syntax: DKBBFTH - energetically coupled disc - Comptonised corona model. Specific notes: "External absorption 
column in units of 10^^ atoms cm^, 'un-Comptonised disc temperature (keV), '^external radius of corona, ''inner radius of the accretion disc, ''photon index, 
^temperature of the Comptonising corona, ^ t is not a fit parameter of the model, but is derived from F and kTg,. 



its innermost stable orbit. This assumption requires that the opti- 
cally thick corona does not intercept our line of sight to the inner 
disc and that the underlying disc spectrum is independent of that 
of the corona (see Kubota & Done 2004). Yet both assumptions 
are probably flawed. An optically thick corona could easily mask 
the innermost regions (depending on its geometry: Kubota & Done 
2004). Secondly, we must consider the energetics of both compo- 
nents. Both the disc and the corona must be ultimately powered by 
gravitational energy release. If we are observing a powerful corona, 
this implies less energy is available for heating the disc (Svensson 
& Zdziarski 1994). The only code currently available that enables 
us to explore the effect of relaxing these assumptions is DKBBFTH 
(Done & Kubota 2006). This was designed to model the extreme 
very high state spectra seen in BHBs such as XTE J1550-564 and 
GRS 1915+105. It incorporates the energetic disc-corona coupling 
model of Svensson & Zdziarski (1994), which assumes that the 
corona extends over the inner disc from Rin to Rt, taking a frac- 
tion / of the gravitational energy available at these radii. Only the 
remaining fraction (1 — /) is available to power the inner disc emis- 
sion, so the inner disc is cooler, but more importantly, less luminous 
than it would be if the corona were not present. This distorted inner 
disc spectrum is the source of seed photons for the Comptonisation 
(Done & Kubota 2006). 

This model has five free parameters (only one more than in 
the purely phenomenological MCD plus power law model and the 
same as for all the disc plus Comptonisation models). The first two 
are the temperature of the Comptonising plasma, kTe, and its opti- 
cal depth (parameterised by the asymptotic spectral index, F). The 
disc is assumed to extend from some inner radius, Ri„, which is 
given by the overall normalisation of the model, and with a tem- 
perature distribution in the limit of no corona being present of 
T(i?)o = Tin{RIRin)~^'^ i.e. as in the diskbb model. How- 
ever, the model assumes that the corona extends homogeneously 
(constant temperature and optical depth) in a slab over the disc, 
from its inner radius to Rt (in units of the Schwartzschild radius, 
Rs = 2GM/c?). For all radii between this and Rin, the calcu- 
lated disc temperature is T[R) = r(i?)o(l - ff^'^ where / is 
the fraction of power dissipated in the corona (assumed constant 
with radius) which is self consistently calculated via iteration from 



the coronal spectral parameters (Done &l Kubota 2006). The result- 
ing disc luminosity underneath the corona is reduced by a factor 
(1 — /), but only a fraction e^^ of this is seen directly, with the re- 
mainder being Compton scattered by the corona (Done & Kubota 
2006). 

Despite being more physically constrained, the model gives 
an equivalently good fit to the data as the disc plus Comptonisation 
models discussed above. Figure|7]shows this fit for Ho IX X-1, with 
the outer (un-Comptonised) disc emission (red solid line) dominat- 
ing at soft energies. The inner disc emission (red dotted line) is 
much lower than in standard disc models as much of the energy is 
powering the corona. These form the seed photons for the Compton 
scattering, but since this scattering is in an optically thick corona, 
these photons are not seen as almost all of them are upscattered to 
form the Comptonised spectrum (solid green line). 

Figure [8] shows the model fit to all our sources, with the pn 
data corrected for absorption to show the intrinsic spectral shape. 
This clearly shows the difficulties in interpreting the parameters 
from simple, phenomenological spectral fits. Again, using Ho IX 
X-1 as an example, the data have an inflection at soft energies that 
characterises the soft excess. The characteristic energy of this fea- 
ture is interpreted as the peak of the disc emission in the MCD plus 
power-law model (giving a very low temperature: hence IMBHs), 
while the high energy break forms the peak of the disc temperature 
in single component disc models (giving a very high temperature: 
too extreme even for stellar mass black holes). In these coupled 
disc-corona models, the 'true' disc temperature is not given by ei- 
ther of these observed features! The high energy break is from the 
very low plasma temperature of the Comptonising region, while 
the low energy inflection occurs when the outer, un-Comptonised 
disc emission starts to dominate the spectrum. The inferred intrin- 
sic inner disc temperatures, recovered from the assumption of disc- 
corona energy partition, range from 0.3 < fcTdisc < 1.2 keV with 
8/12 giving fits where fcTdisc > 0.5 keV, which lies in the stellar 
mass black hole regime. This can be seen explicitly by converting 
the inferred inner radius into black hole mass assuming that the 
disc extends down to the last stable orbit around a Schwarzschild 
black hole at Ri^ = 6GM/c^ = 8.9M. All these 'hot' ULXs 
(fcTdisc > 0.5 keV) have inferred black hole masses < 100 Mq, 
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consistent with stellar mass black holes, and consequent derived 

Lx/LEdd ^ 1- 

However, a third of these source spectra are best fit by fcTdisc 
< 0.5 keV, hence giving black hole masses in the range 80-430 Mq 
and 5MZ;-Eddington accretion rates assuming that the disc extends 
down to the last stable orbit. We propose instead that this very cool 
temperature and hence large radius is not associated with the direct 
disc emission but instead arises as the accretion rate increases be- 
yond that of standard super-Eddington accretion. In a super-critical 
accretion regime, the inner regions of the disc and corona may be 
blown off in the form of a wind, forming an optically thick photo- 
sphere out to large radii (e.g. Poutanen et al. 2007). We return to 
this point in the next section. 



5 DISCUSSION: ULTRALUMINOUS X-RAY SOURCES, 
SUPER-EDDINGTON ACCRETION AND THE 
ULTRALUMINOUS STATE 

We have examined twelve of the highest quality ULX data sets cur- 
rently available in the public archives of the XMM-Newton tele- 
scope. Initial characterisation of the spectrum of these sources re- 
veals the presence of both a soft excess and a break at higher ener- 
gies (within the XMM-Newton band pass). The existence of curva- 
ture at lower energies is apparent in many of the standard BHB 
accretion states, and is regularly fit by an accretion disc (Done, 
Gierlinski & Kubota 2007). It is this standard practice which led 
to the suggestion that these objects were intermediate mass black 
holes, due to the low temperature of the apparent disc emission. 

However, the detection of a break or curvature at higher en- 
ergies (~ 5 keV) brings new insights. Such a feature has been 
observed in ULXs previously (e.g. Foschini et al. 2004; Feng & 
Kaaret 2005; SRW06; Miyawaki et al 2009), but here we show 
that it is nearly ubiquitous in the highest quality spectral data, with 
> 10,000 counts. No such break at these low energies is observed 
in the high energy tail of any of the standard accretion states ob- 
served in BHBs (Remillard & McClintock 2006; Done, Gierlinski 

6 Kubota 2007f| so this challenges the basic assumption of the 
IMBH model, which is that we are observing standard accretion 
states that are scaled with the mass of the compact object. Thus we 
must consider the alternative, that we are observing a different ac- 
cretion state than those generally seen in BHBs. A new accretion 
state would require us to be observing a different mass accretion 
rate with respect to Eddington than those seen in the standard states. 
Since these span (~ 10"'^ - 1) LEdd (McClintock & Remillard 
2006) then this is most likely a super-Eddington state, so requires a 
stellar mass rather than IMBH accretor. 

We therefore suggest that a new observational state should be 
defined based on the characteristic signatures of ULXs. The ultra- 
luminous state is one in which we observe a new combination of 
observational signatures; both a cool disc and a break or roll over at 
high energies in the band pass of the XMM-Newton telescope. We 
caution that high quality data is required for a clear identification; 
the high energy break in particular is difficult to identify with less 
than 10, 000 counts in the X-ray spectrum. 

In order to explore the physical origins of this state we applied 
more physically motivated models to the data. One theory that has 
been presented to explain such an accretion state is the 'slim disc' 

^ Whilst there is a break in the high energy tail of the low/hard and very 
high states, it is at energies > 100 keV and 20-30 keV, respectively cf. 
McClintock & Remillard (2006). 



model, where advection of radiation suppresses the emitted disc 
luminosity of the innermost regions (Abramowicz et al. 1988). We 
apply a simplifed model often used for such spectra, where the tem- 
perature profile in the disc is assumed to be T oc r~^, with p a free 
parameter rather than fixed at 0.75 as for standard discs (Watarai 
et al. 2001). Our data give 0.4 < p ^ 0.6, similar to the p = 0.5 
expected for advection dominated discs, but the derived inner disc 
temperatures are unrealistically high (Tin ^ 6 keV) for one third of 
our sample. Figure |4] illustrates this for NGC 1313 X-1, where T 
~ 8 keV. Plainly, the disc temperature is set by the presence of the 
high energy break, but such high inner disc temperatures are not 
expected even for stellar mass black holes (Ebisawa et al 2003). 
However, we caution against interpreting these parameters physi- 
cally as the model fails to account for the inflection that occurs at 
~ 2 keV. These slim disc models cannot simultaneously produce 
both the soft excess at low energies and the high energy break seen 
in the high quality data used here, though this deficiency may be 
hidden in lower signal-to-noise data. 

In order to explore the nature of these sources in more detail, 
we consider a combination of a disc plus Comptonisation due to 
its successes in describing the emission of other accretion-powered 
BHB systems. Two different Comptonisation models are applied 
to the data, COMPTT and EQPAIR, giving two slightly different 
approximations to thermal Compton up-scattering of accretion disc 
photons (see Figure |6]l. Irrespective of which model is used, we 
find evidence for a cool, optically thick corona, where the high en- 
ergy break is set by the electron temperature of this Comptonising 
plasma. The parameters of this corona are rather different than any- 
thing observed in any of the standard accretion states of BHBs. It 
has lower temperature and higher optical depth than even the most 
extreme optically thick corona seen in the very high state (Done & 
Kubota 2006). Such material must block our view of the inner disc 
(Kubota & Done 2004), and most likely also changes the energy 
dissipated in the disc (Done & Kubota 2006). 

We attempt to recover the intrinsic disc spectrum by mod- 
elling a corona over the inner disc that both Comptonises the in- 
ner regions of the disc and is energetically coupled to it. Both 
corona and disc are ultimately powered by gravity, so increasing 
the power dissipated in the corona at a given radius must mean that 
the disc is less luminous than expected at that radius (Done & Kub- 
ota 2006; Svensson & Zdziarski 1994). Again the parameters indi- 
cate a more extreme version of the very high state, but unlike the 
phenomenological disc plus Comptonisation models, these coupled 
disc-corona models allow us to infer what the source would have 
looked like without any corona. 

The spectral energy distributions derived on the basis of this 
model can be put into a potential sequence of ULX spectra, as 
shown in Figure [8] The first class are those which have spectra 
which increase monotonically in vfv, with maximum power out- 
put at the energy given by the high energy break. All these sources 
give a 'hot disc' (fcTin > 1 keV) in the canonical MCD plus power- 
law fits, with the 'power-law' producing the additional flux at the 
softest energies. This is probably due to the MCD model provid- 
ing a poor description of the broader spectra expected from more 
realistic disc models (Done & Davis 2008; Hui & Krolik 2008). 
However, these spectra still look fairly similar to the standard disc 
spectra seen in the disc-dominated state (NGC 2403 X-1; M81 X- 
6 and M33 X-8). Small amounts of emission from a hot corona 
can also contribute to the spectrum at the highest energies (see e.g. 
the spectral decompositions for the most luminous states on XTE 
J1817-330 in Figure 4 of Gierlinski, Done & Page 2009). We note 
that all of these ULXs have luminosities ^ 3 x 10''^ erg s^^, so 
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Figure 8. XMM-Newton EPIC pn data (black) for all sources in our sample, absorption-corrected and deconvolved with DKBFTH (shown in blue). The 'tnte' 
disc spectrum is over-plotted in red, this is the disc spectrum that would be observed in each case if the corona was removed. It is evident from these spectral 
plots that we are observing a variety of spectral shapes. The first four objects (NGC 2403 X-1 , M81 X-6, M33 X-8 and NGC 1313 X-2) appear very disc-like 
in structure and could be representative of the high or very high state. As we move further down the plots an inflection begins to emerge at ~ 2 keV, signifying 
a break from the standard sub-Eddington accretion states, which we suggest represents a transition to a new super-Eddington accretion state. As the apparent 
disc temperature cools, the spectrum tips, indicating the possible presence of strong winds enveloping the inner regions of the accretion disc, leading to the 
most extreme cases (Ho II X-1, NGC 5408 X-1 & NGC 55 ULX). 



can be close to Eddington for moderately massive (30 — 5Oi\f0) 
stellar remnant black holes, similar to that found in IC 10 X-1 (see 
introduction). The next source in the sequence, NGC 1313 X-2, has 
a spectrum where the high energy emission seems stronger than ex- 
pected from a disc dominated state, so this could instead be a type 
of very high state, again with a moderately massive stellar remnant 



black hole. Hence, it appears as though the low luminosity end of 
the ULX population could potentially overlap with sub-Eddington 
processes seen in the BHB, albeit for larger black holes. We will 
explore this in a future paper by characterising the properties of the 
Galactic BHBs at high Eddington fractions in the XMM-Newton 
band pass. 
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The next category are those where there is clearly a soft in- 
flection as well as a high energy break, but where the total power 
still peaks at the high energy break (Ho IX X-1; IC 342 X-1; NGC 
1313 X-1). These are the ones where the simple disc plus (broken) 
power-law fits give a cool disc together with a high energy rollover, 
which we identify with a new ultraluminous state. For these data 
the coupled disc-corona models give an intrinsic disc temperature 
that is not given by either of these observed characteristic energies, 
and where the effect of Comptonisation in an optically thick, low 
temperature corona, is most marked. These spectra do not corre- 
spond to any of the known states, but can form from a more extreme 
(higher optical depth, lower temperature) version of the very high 
state corona. These are typically brighter than those in the previous 
class (though there is also substantial overlap in luminosity) so are 
most likely super-Eddington accretion flows. 

There is then a clear observational sequence of spectral 
shapes, through the sources where the ratio of power between 
Comptonisation (the high energy peak) and the outer disc (low en- 
ergy peak) steadily decreases, from NGC 5204 X-1 and NGC 4559 
X-1, to the most extreme systems, namely Ho II X-1, NGC 5408 X- 
1 and NGC 55 ULX. These are the ones where the inferred intrinsic 
disc temperature from the coupled disc-corona model is < 0.4 keV, 
far lower than expected from stellar remnant accretion. 

To understand these, we look first at what happens physically 
to the flow as it approaches and then exceeds the Eddington limit. 
It has long been known that the Eddington limit for a disc is some- 
what different than that for spherical accretion (Shakura & Sun- 
yaev 1973). A thin disc, which radiates at the Eddington limit at all 
radii, has an integrated luminosity L '--^ (1 + in(m))I/Edd where 
m = M/M-Edd is the mass accretion rate scaled to that which 
gives the Eddington luminosity for spherical accretion. However, 
this can only be achieved if excess energy over and above that ex- 
pected from a constant mass inflow rate, radiating with constant ef- 
ficiency, is somehow lost from the system. There are two ways to do 
this, either by changing the mass accretion rate as a function of ra- 
dius through expelling the excess mass via winds (Shakura & Sun- 
yaev 1973; Begelman et al. 2006) or by changing the radiative effi- 
ciency by advecting the photons along with the flow (slim discs, as 
above). Importantly, both can operate simultaneously (Poutanen et 
al 2007), as indeed is shown in the most recent 2D radiation hydro- 
dynamic simulations of super-Eddington accretion flows (Ohsuga 
2006; 2007; 2009; Kawashima et al. 2009; Takeuchi , Mineshige, 
& Ohsuga 2009). 

This then gives a possible framework to interpret our results. 
It is clear that the objects with ~ 1 keV disc emission are proba- 
bly just more extreme versions of the brightest high and very high 
spectral states known from BHBs. The super-Eddington, ultralu- 
minous state sources are then distinguished by their cool, optically 
thick coronae and apparently cooler discs. One obvious source of 
higher optical depth in the corona is the increasing importance of 
winds as the source starts to accrete past the Eddington limit, and 
this mass loading of the coronal particle acceleration mechanism 
leads to lower temperatures of the Comptonising electrons. 

These winds will become increasingly important as the flows 
become increasingly super-Eddington, completely enveloping the 
inner regions of the disc-corona out to an increasing photospheric 
radius (as in SS 433, Poutanen et al. 2007). Hence the ob- 
served temperature decreases in line with the Stefan-Boltzmann 
law (Shakura & Sunyaev 1973; Begelman et al. 2006; Poutanen et 
al. 2007). The outflow is inherently (at least) two dimensional, so 
viewing angle will change the apparent system luminosity (Ohsuga 
2006; 2007). We suggest the most extreme objects seen in our cur- 



rent sample, NGC 4559 X-1, Ho II X-1, NGC 5408 X-1 and NGC 
55 UL^J, are dominated by reprocessing in a wind, and that much 
of their luminosity output is channeled into kinetic energy. 

Thus it seems most likely that we are seeing the ULX tran- 
sit between the brightest high and very high states, through to a 
super-Eddington ultraluminous state which is similar to the very 
high state but with lower temperature and higher optical depth in 
the corona, to a completely new (hyper-accreting) state where the 
wind dominates the spectrum. In none of these states do we require 
the presence of IMBHs to explain the X-ray spectrum; all can be 
explained by stellar mass black holes at high accretion rates, albeit 
perhaps black holes up to a few times larger than those known in 
our own Galaxy. 



6 CONCLUSIONS 

The highest quality data has been collated and utilised to both char- 
acterise the spectra of ultraluminous X-ray sources and to constrain 
their nature. These show that while some ULX (typically the low- 
est luminosity ones) have spectra which are probably similar to the 
high and (especially) the very high state in Galactic BHBs, the ma- 
jority show more complex curvature which can be modelled by a 
cool disc component together with a power-law which breaks/rolls- 
over above ~ 3 keV. This combination of spectral features is not 
commonly present in any of the known (sub-Eddington) Galactic 
BHB states, and we therefore propose these features as observa- 
tional criteria for a new ultraluminous state and identify it with 
super-Eddington accretion flows. 

More physical models for these spectra show they are not well 
fit by (approximate, 'p-free') slim disc models as these cannot si- 
multaneously produce both the soft excess and high energy break. 
Instead, disc plus Comptonisation models give a much better de- 
scription of this complex curvature, indicating that the break above 
~ 3 keV comes from a cool, optically thick corona. This suggests a 
more extreme version of the coronae seen in the very high state of 
BHBs, as might be expected for super-Eddington flows. However, 
such coronae obscure the inner disc and alter its energetics (Done 
& Kubota 2006), so we model these effects to recover the intrin- 
sic disc temperatures. Many of these are in the range expected for 
stellar remnant black holes, showing that the apparent cool disc 
temperature derived from simple disc models does not require an 
IMBH. However, there are some objects where the recovered disc 
temperature (corrected for the corona) is cooler than expected for a 
stellar remnant black hole. We suggest these are most likely to rep- 
resent the most extreme super-Eddington accretion flows, where the 
wind from the accretion disc becomes so powerful that it envelops 
the inner disc out to a large photospheric radius, producing the cool 
spectral component. 

There are occasional spectra from BHBs at the highest lumi- 
nosities which are indeed better described by optically-thick Comp- 
tonisation (GRO J1655-40 and GRS 1915+105: Makishima et al. 
2000; Middleton et al. 2006; Ueda et al. 2009). Similarly, the high- 
est Eddington fraction AGN such as RE J 1034+396 (Middleton et 
al. 2009) and RX JO 136-35 (Jin et al. 2009) also show such spectra. 
We suggest that all these sources are in this new super-Eddington 
ultraluminous accretion state. It now appears that ULXs, rather than 

^ The relatively low luminosity of NGC 55 ULX may be an effect of its 
disc being close to edge-on to our line-of-sight, evidence for which comes 
from the dipping behaviour seen in its X-ray light curve (Stobbart et al. 
2004). 
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revealing a new population of IMBHs, are providing us with a tem- 
plate for accretion at super-Eddington rates. This will have wide 
applications across many areas of astrophysics, ranging from stel- 
lar formation to the growth of QSOs. Further studies of ULXs that 
provide a deeper understanding of this new and crucially important 
accretion regime are therefore imperative. 
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